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A generalised model for interacting self-avoiding trails on the square lattice is presented and 
studied using numerical transfer matrix methods. The model differentiates between on-site double 
visits corresponding to collisions, and crossings. Rigidity is also included in the model. The model 
includes the Nienhuis O(n=0) model and the interacting self-avoiding trail model as special cases. 
It is shown that the generic type of collapse found is the same as the pure interacting self-avoiding 
trail model. 
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Lattice walk models have been used for many decades 
as coarse-grained models of polymers in solution[l|. The 
quality of the solvent is modelled by introducing short- 
ranged interactions 0, Whilst the interacting self- 
avoiding walk, with an attractive interaction associated 
with nearest-neighbour non-consecutively visited sites, is 
the most studied lattice walk model, models in which an 
attractive interaction is associated with multiply visited 
sites have also attracted attention. These models fall into 
two types: the Vertex-Interacting Self- Avoiding Walk 
(VIS AW) model [1] and the Interacting Self- Avoiding 
Trails model [|. The VIS AW is a bond-avoiding walk 
which is allowed to revisit sites, but is not allowed to in- 
tersect itself, whereas the ISAT is allowed to intersect at 
a lattice site. 

The VISAW is rather well understood, partly because 
it corresponds to the n=0 limit of the 0{n) model in- 
troduced by Blote and Nienhuis [4j. The ISAT model, 
on the other hand, has given rise to many apparently 
contradictory results;^. In recent work we presented ev- 
idence that the collapse transition in the ISAT trail has 
the same correlation length exponent, v, as the VISAW 
but a different entropic exponent 7[3]- The apparent con- 
tradiction with the results of Owczarek and Prellberg's'l 
is resolved by the realisation that the relation v — l/df 
between the v exponent and the fractal dimension does 
not always holdQ. 

An interesting question arises as to the generic nature 
of the ISAT point. If the parameter space is enlarged to 
enable one to extend the collapse point to a line, is the 
ISAT the generic type of collapse, or does one recover 
a more conventional 0-type collapse? This question was 
recently examined on the triangular latticetld], where the 
construction of the lattice leads naturally to two different 
attractive interactions; on the triangular lattice one may 
revisit a site once or twice without visiting the same bond 
twice. In this study Doukas et al found that the 
ISAT collapse corresponded to a special point separating 
an isotropic collapse phase, a crystalline phase and the 
swollen phase. In this phase space the ISAT collapse 
extended into a line of ^-points. 

In this brief report we study a generalised interacting 
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FIG. 1: The generalised interacting self-avoiding trail model 
showing the vertex crossings, weighted with a Boltzmann fac- 
tor Tx, and the vertex collisions, weighted with a Boltzmann 
factor Tc. A weight p is associated with straight segments 
and a fugacity K controls the length of the walk. The walk 
is shown on a strip of width L = 5 and length Lx = 10. 



self-avoiding trail model on the square lattice to exam- 
ine how the collapse varies as the model is moved from 
the standard interacting trail model. In our model this 
may be achieved either by differentiating between cross- 
ings and collisions, or by introducing a rigidity (weighting 
straight sections with respect to corners). 

The model studied here is defined as follows: consider 
all random walks on the square lattice which do not visit 
any lattice bond more than once. Doubly visited sites 
may correspond to either crossings or "collisions" , cross- 
ings are assigned an attractive energy —ex^ whilst col- 
lisions are assigned an attractive energy —£c- We also 
introduce a penalty Es for straight sections. The parti- 
tion function for the model is 



Z = 



E 

walks 



T^N 



(1) 



where K is the step fugacity, t^/c = ^'^viP^x/c), P — 
exp(— /3es), N is the length of the walk, and is the 
number of crossings and Nc is the number of collisions. 

This partition function may be calculated exactly on a 
strip of length ^ oo and of finite width L by defining 
a transfer matrix T. If periodic boundary conditions are 
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assumed in both directions, the partition function for the 
strip is given by: 



Zl = lim Tr (T^ 



(2) 



The free energy per lattice site, the density, and correla- 
tion length for the infinite strip may be calculated from 
the eigenvalues of the transfer matrix: 
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(3) 
(4) 
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where Aq and Ai are the largest and second largest (in 
modulus) eigenvalues. 

It is expected that Z, p and ^ should have the following 
scaling forms close to the critical fugacity (for fixed t): 

Z ^ \K~K,\-'\ (6) 
^ ^ \K-K,r, (7) 
Pl{K) = p^(K) + L'/''~^p{\K~K,\L'/n. (8) 

Z corresponds to the high temperature expansion of the 
susceptibility of an equivalent magnetic model, hence the 
use of the exponent 7. 

These scaling properties enable estimates of the criti- 
cal lines to be calculated using a phenomenological renor- 
malisation group method. For example a critical point 
estimate for a pair of lattice widths L and U is given by 
the solution of the eg nation [l3|: 



L 



L' 



with estimates of the critical exponent v given by: 



log (^/^ 
log(i/i') 



- 1. 



(9) 



(10) 



For a more detailed discussion of the transfer matrix 
method, the reader is referred to the article of Blote and 
Nienhuis [I]- 

We first study the flexible generalised ISAT model. By 
flexible we mean that p — 1 and the straight sections 
and corners are equally weighted. The phase diagram 
was calculated in the plane where the average length of 
the walk diverges, which occurs when Ai — Aq, or equiv- 
alently when ^ ^ 00. This (critical) surface may then 
be estimated in two different ways - either by setting 
Ai = 1 (since Aq = 1), or by using the Phenomenological 
RG equation In both cases the collapse transition 
may be found in a variety of ways, but commonly one 
looks at the variation of v as Tc (or Tx) is varied for fixed 
X = Tx/tc- The crossings of v give an estimate of the 
collapse transition. One may calculate v either from the 
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FIG. 2: (Color online) Phase diagram in the x, Tc plane for 
the fully flexible generalised ISAT model {p = 1). The point 
marked A corresponds to the pure VISAW model and the 
point marked B corresponds to the pure ISAT model. 
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FIG. 3: (Color online) The v exponent as a function of Tc 
calculated for the flexible generalised ISAT model using the 
scaling of p setting Ai = Aq — 1. The dotted line corresponds 
tou^ 12/23. 



scaling behaviour of p using ([5]), useful in the case we 
estimated the critical surface by setting Aq = Ai = 1, 
or by the Phenomenological RG equation pO|) . The two 
methods give superposable results for the phase diagram 
plotted in the critical surface, shown in Figure[5] The es- 
timates of v calculated using the two methods are shown 
in Figures [3] and S] 

Rememberin g th at the VISAW model is expected to 
have v = 12/23 111|, and it is conjectured that the ISAT 
model also has f = 12/23[3|, we may interpret the v re- 
sults as indicating that the entire collapse line to the left 
of the ISAT collapse point (B) is in the ISAT univer- 
sality class, whilst at some stage the collapse transition 
becomes first order. It seems that the ISAT behaviour 
extends beyond a; = 1 {x/{l -I- a;) = 0.5). Both the VI- 
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FIG. 4: (Color online) The v exponent as a function of Tc 
calculated for the flexible generalised ISAT model using phe- 
nomenological RG (equations (O and (|10p ). The dotted line 
corresponds to v — 12/23. 



SAW model and the ISAT model are believed to have 
the same value of ly, whilst the value of 7 was found to 
be different (7 — 53/46 in the first case[il| and conjec- 
tured to be 7 = 22/23 in the second Q)- This difference 
is similar to that seen between the 9 point model on the 
regular lattice where 7 = 8/7[l3| and on the Manhattan 
lattice 7 = 6/7 [l3|- Both the Manhattan lattice 9 point 
model and the ISAT model have the particularity that if 
the walks are grown dynamically, the construction of the 
walk can only fail if the walk meets its starting point. 
This feature is common to the whole collapse line except 
for the VIS AW collapse point (A) {x = 0); since cross- 
ings are forbidden at this point, when walks are grown 
they may find themselves trapped. This suggests that 
the generic behaviour along the collapse line might be 
that of the ISAT, though this is not totally clear from 
the transfer matrix data available. 

On the square lattice, the other obvious means of mov- 
ing away from the ISAT collapse is to include a rigidity 
to the walk. Here we chose to weight the straight sections 
relative to the corners with a weight p. This choice was 
made to coincide with the previous work of Blote and 
Nienhuis[l]- 

The phase diagram in the critical surface is calculated 
as before, and shown in Figure [5] in the Tc,p plane for 
X = 1. The limit p — excludes the possibility of cross- 
ings, since a crossing necessarily requires two straight 
portions. At this point the model may be mapped into 
the interacting self-avoiding walk on the Manhattan lat- 
tice, which, as discussed above, is known to be 0-like 
{1/ = 4/7). Following from the results for the ISAT on 
the triangular lattice, where it was found that the ^-type 
behaviour was more generic, we may expect that the col- 
lapse transition between (A) and the ISAT collapse point 
(B) might fall into a 9-like universality class (0-like since 
the 7 exponent may differ from the 9 point value 8/7). 
The results, whilst not conclusive, do seem to indicate 



FIG. 5: (Color online) Phase diagram in the p, Tc plain for 
the semi-flexible ISAT model {x = 1). The point marked 
A corresponds to the 0{n = 0) Nienhuis-Blote model with 
only corners. This point then maps onto the point on the 
Manhattan lattice or on the L lattice (see text). The point 
marked B corresponds to the pure ISAT model. 
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FIG. 6: (Color online) The 1/ exponent as a function of Tc cal- 
culated using phenomenological RG (equations ([9]) and (|10p ') 



that the value calculated for the ISAT model extends 
into a plateau as the lattice widths increase (see Fig- 
ure[6|), leading to the conclusion that the ISAT behaviour 
is again in this case more generic. Again the results indi- 
cate that the plateau extends beyond the ISAT collapse 
point, before the collapse transition becomes first order 
for some value of p > 1 . 

In this brief report we presented transfer matrix results 
for a generalised model for interacting self-avoiding trails 
in order to investigate the stability of ISAT collapse to 
varying the model parameters. On the triangular lattice 
Doukas et al [l0| presented convincing Monte Carlo re- 
sults which showed that a collapse transition similar to 
the ISAT collapse was a special point in the phase di- 
agram, and that generically the collapse transition was 
either in the 9 universality class or first-order. The first 
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type of collapse separates a swollen walk from a dense 
isotropic phase and the second separates the swollen walk 
from an anisotropic crystalline collapsed phase. On the 
square lattice the ISAT collapse transition is an extended 



line separating the swollen phase from the crystalline 
phase. The collapse transition becomes first order when 
crossings dominate over collisions, which happens either 
for large enough x = Tx/tc or large enough p. 
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